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ABSTRACT 

A principal component analysis of cosmic microwave background (CMB) anisotropy 
measurements is used to investigate degeneracies among cosmological parameters. 
The results show that a degeneracy with tensor modes - the 'tensor degeneracy' - 
dominates uncertainties in estimates of the baryon and cold dark matter densities, 
ojb — ^bh^, <^c = ^ch^^ from an analysis of CMB anisotropies alone. The principal 
component analysis agrees well with a maximum likelihood analysis of the observa- 
tions, identifying the main degeneracy directions and providing an impression of the 
effective dimensionality of the parameter space. 

Key words: cosmic microwave background- cosmology: miscellaneous. 
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1 INTRODUCTION 

Since the discovery of CMB anisotropies by the COBE team 
(Smoot et al. 1992) there has been rapid progress on the ob- 
servational front, culminating with publication earlier this 
year of evidence for multiple acoustic peaks in the CMB 
power spectrum (Netterfield et al. 2001; Lee et al. 2001; 
Halverson et al. 2001). It has long been known that accurate 
measurements of the CMB anisotropies can be used to esti- 
mate parameters characterising the primordial fluctuations, 
the geometry of the Universe and its matter content {e.g. 
Jungman etal 1996; Bond, Efstathiou & Tegmark 1997; Zal- 
darriaga, Spergel & Seljak 1997). In fact, the most recent 
measurements paint a gratifyingly consistent picture com- 
patible with the simplest models of inflation {i.e. a spatially 
flat Universe with scale-invariant adiabatic fluctuations ). 
Furthermore, the derived value of the baryon density cub ap- 
pears to be consistent with the value ui, = 0.020 ± 0.002 in- 
ferred from primordial nucleosynthesis and deuterium abun- 
dance measurements from quasar absorption line spectra 
(Buries, NoUett & Turner 2001, and references therein). 

It has also long been known that there are significant 
degeneracies amongst cosmological parameters estimated 
from CMB anisotropies, i.e. parameter combinations exist 
that produce nearly identical CMB power spectra (Bond 
et al. 1994; Efstathiou & Bond 1999, hereafter EB99). The 
best known is the geometrical degeneracy between the mat- 
ter and vacuum energy densities, flm and Ha, and the cur- 
vature rifc = 1 — flm — f^A. This degeneracy is almost exact 



Here h is Bubble's constant Ho in units of lOOkms ^Mpc 



and precludes reliable estimates of either Qa or the Hubble 
parameter h from measurements of the CMB anisotropies 
alone. The existence of parameter degeneracies means that 
the best fitting parameters and their errors can be extremely 
sensitive to the chosen parameter set {e.g. whether the Uni- 
verse is assumed to be spatially flat) and to adopted 'prior 
distributions' {e.g. observational constraints on the Hub- 
ble constant). This complicates the interpretation of CMB 
anisotropy parameter studies and the intercomparison of 
limits determined by different authors. 

Almost all CMB parameter analyses (with the note- 
able exception of the work of Tegmark and collaborators, 
(Tegmark 1999; Tegmark, Zaldarriaga & Hamilton, 2001; 
Wang, Tegmark & Zaldarriaga, 2001) have ignored a ten- 
sor component (see for example, Lange et al. 2001; Jaffe et 
al. 2001; de Bernardis et al. 2001; Pryke et al. 2001; Stom- 
por et al. 2001). This special case is certainly interesting be- 
cause a wide class of inflationary models predict a negligible 
contribution from tensor modes. Indeed it has been argued 
persuasively that the absence of tensor modes is generic to 
any model in which the inflaton potential is related to the 
Higgs sector of a grand unified theory (Lyth 1997). 

However, so little is known about inflation (if indeed 
inflation occured) that it may be dangerous to neglect a 
tensor mode, particularly if some cosmological parameters 
of interest are sensitive to tensor modes. In fact, in single- 
field inflation models the relative amplitudes of the tensor 
and scalar tensor modes (r) and their spectral indices {nt 
and ria) are related to the infiaton potential and its first 
two derivatives {e.g. see Hofl'man & Turner 2001 for a re- 
cent discussion) . In this class of model the relation between 
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Us and rit is model dependent while the relation between nt 
and r depends only on the validity of the slow-roll approxi- 
mation. However, even the latter relation can be violated in 
multi-field infiationary models. These relations may also be 
violated in some superstring inspired models (see e.g. Lidsey, 
Wands & Copeland 2000, and references therein). Examples 
of the latter include the pre-big bang model of Veneziano and 
collaborators {e.g. Buonnono, Damour &Veneziano, 1999) 
which it is argued can lead to 'blue' scalar spectral indices 
(wa > 1), and the ekpyrotic scenario (Khoury, Ovrut, Stein- 
hardt & Turok 2001) which produces a strongly blue tensor 
mode spectrum. 

The point of view taken in this paper is to define a 
minimal parameter set on the assumption that the primor- 
dial fluctuations are Gaussian, adiabatic and featureless {i.e. 
defined by power-law spectral indices). Thus the minimal 
model is specified by 9 parameters: four parameters specify- 
ing the amplitudes and spectral indices of scalar and tensor 
components (Q, rio, Us and nt, see Section 2 for more pre- 
cise definitions), four parameters defining the matter content 
and curvature of the Universe (lj;,, uJc, ^a, Q.^) and a sin- 
gle parameter Topt quantifying the optical depth to Thom- 
son scattering since the Universe was reionized. However, as 
pointed out by EB99, including a tensor component with no 
assumed constraints between the parameters ris, nt and rio 
introduces an new major degeneracy between cosmological 
parameters that we will call the tensor degeneracy in this 
paper. This degeneracy has a dramatic effect on the permit- 
ted ranges of some parameters, in particular the baryon and 
cold dark matter densities Ut, and uJc. 

One approach to these degeneracies is to apply brute- 
force maximum likelihood analysis to a large parameter 
set (see the papers by Tegmark and collaborators). Here 
we show that the observational data has now improved to 
the point that a simple principal component analysis of the 
Fisher matrix defined by the data identifies the major de- 
generacy directions (geometrical and tensor) and provides a 
useful first approximation to the likelihood function. Using 
the principal components it is easy to analyse the correla- 
tions between physical parameters introduced by the tensor 
degeneracy and to assess the effects of introducing external 
(non-CMB) constraints on the parameters. 



2 PRINCIPAL COMPONENT ANALYSIS 

We use the compilation of band power estimates ATg 
and their covariance matrix Cbb' computed by Wang et 
al. (2001, hereafter WTZOl) from 105 CMB anisotropy mea- 
surements. These band power estimates include a model for 
calibration and beam errors (see WTZOl for further details). 
Each band power estimate is related to the power spectrum 
Ct of the CMB anisotropics by 



ATS = ±0- 



2-K 



< 




10 20 



Figure 1. The points show band-averaged observational esti- 
mates of the CMB power spectrum from WTZOl plotted against 
multipole inde using a log-linear abscissa. The error bars show 
itlcr errors. The line shows the CMB power spectrum for the 
fiducial inflationary model discussed in the text. 



cial model|^ with the following parameters: uji, = 0.020, 
LUc = 0.13, h = 0.7 (fib = 0.04, f7c = 0.26), 9.k = 0, 
Topt ~ 0.1, ris — 1, nt — 0, rio = 0.2. These parameters 
provide an extremely good fit to the observations and are 
very close to the concordance values determined by WTZOl 
and the author from a full likelihood analysis. The addition 
of a small tensor component has little effect on the fit shown 
in Fig. 1, but is introduced to regularize the Fisher matrix 
of equation (^). 

Given the covariance matrix of the observations, we can 
form the Fisher matrix for the (mean subtracted) parameter 
set {si}: 



1 

BB 



dATl; aAr|, 



dsi 



dsj 



(2) 



^^(£+ l)CfW^s(£) 



(1) 



In defining the parameter set st, we use Incob and Incite rather 
than ujb and lJc and iln = fife — 0.286f2A rather than Q,k- 
The latter expresses the geometrical degeneracy (see EB99) , 
since if all parameters other than flk and are held fixed, 
the condition Sflo ~ preserves the positions of the acous- 
tic peaks for small variations of the parameters about those 
of the fiducial model. The normalization parameter Q is de- 
fined following EB99 so that is the mean band power 
(X/f<i500 ^(^+ 1)C'«)) of model relative to that of the fidu- 
cial model plotted in Fig. 1. With this definition, Q is ob- 
servationally well constrained and is largely decoupled from 
variations in cosmology and the optical depth Topt (unlike 
measures of the amplitude related to low multipoles). The 
tensor to scalar ratio rio is defined so that Cio = rioCfo. 
With these definitions, we can compute the Fisher matrix 
(□) from derivatives of the power spectra Cf and Cj, as 



where Wb is the window function for each band power 
(also computed by WTZOl). These band-power estimates 
are plotted in Fig. 1. Also plotted in this figure is a fidu- 



t The CMB power spectra in this paper were computed using the 
CMBFAST code of Seljak & Zaldarriaga (1996). 
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Figure 2. Various projections of tlie likeliliood functions chosen to illustrate the geometrical and tensor degeneracy. One, two and 
three-sigma likelihood contours are plotted in each panel. The panels to the left were computed from the Fisher matrix (equation ^ and 
the panels to the right show 'pseudo-marginalized' contours computed from a full likelihood analysis of the data described by Efstathiou 
et al. (2001). The solid lines in the upper panel show the approximate locus of the geometrical degeneracy computed from the constraint 
Qo = constant and from the constraint equation (^. The dashed lines in the upper and lower panels show the degeneracy directions 
defined by the component Jfg and Xg respectively. The filled circles show the parameters of the nearly degenerate models plotted in Fig. 
3. The crosses in figures 2b show the positions of the peaks in the likelihood function. 



in standard analyses of parameter forecasting (Jungman et 
al. 1996, Bond et al. 1997). Alternatively, we could compute 
the Hessian matrix by estimating second derivatives of the 
likelihood function around the maximum likelihood value ei- 
ther directly or by summing over first and second derivatives 
of ATg. In most situationa using either the Fisher and Hes- 
sian matrixes should give closely similar results, and in fact, 
the second derivatives of ATg are often ignored in order to 
stabilise numerical evaluations of the Hessian (see e.g. Press 
et al. 1992). 

Having computed the Fisher matrix, we diagonalise it, 

F = UAU'^, A = diaff(Ai,A2,...A]v). (3) 

The matrix U defines a set of principal components, X, i.e. 
orthogonal linear combinations of the original parameters, 

X = U^s, (4) 

such that the variance of the component Xi is equal to 1/Ai. 
(See BE99 for an application of principal component analysis 
to the MAP and Planck satellites^. 

The eigenvalues and components of for our chosen 
set of variables and fiducial model are listed in Table 1. The 
principal components have been ordered by their expected 
variance so that Xi is the best determined parameter and Xg 

■t Descriptions of these satellites can be found on the following 
web pages: http://astro.estec.esa.nl/SA-general/Projects/Planck 
and http://map.gsfc.nasa.gov 



is the worst. The last line in the table lists the predicted vari- 
ances of the parameters Si, {yj {F~l^)). As can be seen from 
Table 1, the present data can be used to constrain three pa- 
rameter combinations well and three extremely poorly, with 
the remaining constrained at intermediate levels of accuracy. 
The best and worst determined principal components have 
a straightforward interpretation. X\ has a high weight from 
Q.D, and so, in effect, measures the positions of the acoustic 
peaks; X2 has a high weight from the overall amplitude of 
the spectrum Q\ X3 has high weights from the scalar spec- 
tral index Us and Topt and provides a measure of the shape of 
the fiuctuation spectrum. The component Xg accounts for 
almost all of the variance of f^A, rit and rio, and describes the 
geometrical degeneracy and the extremely poor constraints 
that the present data place on tensor modes. The compo- 
nent Xs describes the tensor degeneracy and accounts for 
the almost all of the variance of the parameters Inoit, \nujc, 
Us and Qd. Together, the components Xq and Xg accounts 
for almost all of the variance of all of the parameters with 
the exception of the amplitude Q. 

The principal component analysis therefore suggests 
that acceptable models lie on a plane within the 9 dimen- 
sional parameter space. This is illustrated in Fig. 2a which 
shows likelihood contours in projections chosen to illustrate 
the geometrical and tensor degeneracies. The upper pan- 
els show the geometrical degeneracy in the parameter pairs 
Ha — fife and h — Qa. The dashed lines show the degeneracy 
directions defined by the component Xq. The solid line in 
the JIa — f^fc plane shows the constraint SQd = 0. The solid 
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Table 1: Principal components for data of Fig. 1 
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line in the h — Oa plane is computed from the constraint 
equation 

with LOb, uJc and Qd set to the values of the target model. 
The dotted lines in the lower two panels show the degeneracy 
directions defined by the component Xg. 

The Fisher matrix and the associated principal compo- 
nent analysis is approximate and it is not obvious a priori 
how well they decribe degeneracies especially for parame- 
ter values that are quite a long way from those of the tar- 
get model. Fig. 2b shows the analogous contours computed 
from a full likelihood analysis of the data of Fig. 1 (see the 
companion paper of Efstathiou et al. 2001 for a detailed dis- 
cussion). The agreement is surprisingly good. The general 
directions of the geometrical and tensor degeneracies follow 
those of the Fisher matrix analysis. There are some differ- 
ences however. The low order CMB multipoles offer some 
discrimination of models with high values of f^A via the in- 
tegrated Sachs- Wolfe effect (see EB99), thus models with 
Oa ^ 0.88 arc excluded by the data. The maximum likeli- 
hood analysis shows that the tensor degeneracy allows high 
values of cot but that the likelihood function falls sharply for 
uJb ~ 0.018 (interestingly only just below the value favoured 
from primordial nucleosynthesis deuterium constraints). 

The three filled circles in Fig. 2 show the parameter 
values for the target model and for two nearly degenerate 
models with extreme values of ujb chosen to lie along the 
tensor degeneracy direction. (The device used here, of vary- 
ing low order principal components to produce degenerate 
pairs of models, was used by EB99 to investigate parame- 
ter degeneracies for the MAP and Planck satellites.) One 
model has uJb = 0.042 and the other has = 0.012. The 
models were intentionally chosen to have similar values of 
Qa, consequently the high baryon density model has a low 
value of h which one might argue is incompatible with inde- 
pendent measurements of the Hubble constant (Freedman 
et al. 2001). However, our purpose here is to show how the 
constaints on uJb derived from CMB anisotropics alone are 
weakened when tensor modes are included. The high and 
low baryon density models are compared to the data and to 
the fiducial model in Fig. 3. 

Despite the very different parameter values, the models 
produce almost identical CMB power spectra by construc- 
tion. The exact likelihood analysis (Fig. 2b) shows that the 
model with the high baryon density is compatible with the 
data at about the 2cr level. The low baryon density model 



is formally excluded by the data at a high level of signif- 
icance (> 3a") because it fails to match the height of the 
first acoustic peak. This is a characteristic feature of mod- 
els with a low baryon density. Nevertheless, the diagram is 
interesting because it shows that that the lower limits on 
LUb arc extremely sensitive to any residual systematic errors 
that might affect the peak heights (see e.g. de Bernardis et 
al. 2001, figure 1). 



3 DISCUSSION 

The reader might question the usefulness of the results pre- 
sented in the previous section. Firstly, the principal compo- 
nent analysis provides only an approximate description of 
the parameter degeneracies, whereas they emerge precisely 
from a brute force maximum likelihood analysis. Secondly, 
models at the extreme ends of the ranges allowed by the ten- 
sor degeneracy have unusual parameters (for example, the 
model with uJb = 0.042 in Fig. 3b has a low Hubble con- 
stant and a high value of the scalar spectral index) and so 
are surely excluded by other observational constraints. We 
discuss each of these points in turn. 

(i) Effective dimensionality: The main use of the principal 
component analysis is to assess the effective dimensionality 
of the space of acceptable models within the multidimen- 
sional space defined by the physical parameters s, and to 
assess whether this effective dimensionality is sensitive to 
changes of the parameter set. For example, let us assume 
that we are interested in the values of the parameters ujt 
and LOc. The principal component analysis tells us that most 
of the variance of these parameters is contributed by only 
two poorly constrained principal components. The values of 
these parameters are therefore affected by major parameter 
degeneracies which can only be removed by imposing exter- 
nal constraints or by performing a fundamentally different 
type of CMB exeperiment. For example, the effects of the 
tensor degeneracy on uJb and ujc can be broken by extending 
the CMB measurents to much higher multipoles (sec EB99) 
and/or by setting limits on a tensor component from an 
analysis of a B-type polarization pattern in the CMB (see 
e.g. Kamionkowski & Jaffe, 2000, and references therein). 

(ii) Complementary information: CMB parameter degenera- 
cies can be broken by invoking complementary information. 
A well known example is the combination of Type la super- 
novae measurements with the CMB to breaJf the geometri- 
cal degeneracy (see e.g. de Bernardis et al. 2001 for a recent 
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Figure 3. The data and fiducial model (solid line) of Fig. 1. The dashed lines show CMB power spectra for nearly degenerate models 
with a high baryon density (Fig. 3a) and a low baryon density (Fig. 3b) chosen to lie along the direction of the tensor degeneracy. The 
dotted line in Fig. 3a shows the contribution of the tensor component. The low model with low baryon density in Fig. 3b has a neglible 
contribution from tensor modes. 



analysis). WTZOl break the tensor degeneracy by combin- 
ing the CMB data with a number of constraints including 
estimates of the power spectrum on small scales from obser- 
vations of the Lya forest (Croft et al. 2001) and limits on the 
Hubble constant from HST Hubble Key Project (Freedman 
et al. 2001). The problem here is that it becomes progres- 
sively more difficult to assess whether parameter values are 
affected by systematic errors as more external constraints 
are applied, particularly if the external constraints involve 
complex observations and assumptions. (For example, how 
can we test empirically whether density fluctuations in the 
inter-galactic medium as traced by the Lya forest match 
those of the dark matter?). Even if the parameter values 
are shown to remain consistent as external constraints are 
applied, there is no guarantee that the final combined like- 
lihood distribution will be accurate. Ideally, we would like 
to apply as small a number of external constraints as pos- 
sible using data sets with well controlled errors. The Fisher 
matrix analysis of the previous section offers a guide as to 
which external constraints will be most effective at breaking 
degeneracies. The effects of external constraints can easily 
be assessed using the Fisher matrix: let Pij be the covari- 
ance matrix of the parameters Si from external constraints, 
then the covariance matrix after combining with the CMB 
is 

Cij = {Fij + Pij^) ^■ 



As meiisLiuuyd iu QyLliuu 2, thy models pluLLed in Fig. 3 



were ch boon by oonatruotion to have aimilar valuoa of fljv, do 
it is clear that combining the CMB results with constraints 



on JIa from I'ype la supernovae observations will have little 
effect on the parameters lui, and ujc- Constramts on the Hub- 
ble constant will be more effective in narrowing the range 
of allowed values of cub and ujc- However, from Fig. 2 we can 



see that the tensor degeneracy exhibits a strong correlation 
between ujt, loc, and the scalar spectral index Us- Even a 
relatively poor constraint on the scalar spectral index will 
lead to a sharp tightening of the uji, and luc error contours. 
(A constraint Aus on the spectral index will narrow the er- 
ror on the baryon density to Alncjb ~ (il7ig/!738)Ans, i.e. 
Alnoji, « 1.04Ans .) This suggests that the recent estimates 
of the galaxy power spectrum from the 2dF Galaxy Redshift 
Survey (Percival et al. 2001) will be highly effective at break- 
ing the tensor degeneracy, leading to tight constraints on the 
matter content of the Universe. This is borne out by a de- 
tailed analysis presented in a companion paper (Efstathiou 
et al. 2001). 
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